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i ' Abstract 

We have calculated inclusive one- and two-jet production in photon-photon collisions 
superimposing direct, single resolved and double resolved cross sections for center of mass 
energies of the LEP1, LEP2 and NLC range. The direct and single resolved cross sections 
^ are calculated up to next-to-leading order. The double resolved two-jet cross section 

is calculated only in LO with a k factor estimated from the NLO one-jet cross section. 
Various differential cross sections as functions of transverse momenta and rapidities of the 
jets are evaluated. 
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1 Introduction 



In a recent paper we presented theoretical results for the production of high-pp jets in almost- 
real photon-photon collisions |I[ and compared them with recent experimental data for inclusive 
one- and two-jet cross sections from the two TRISTAN collaborations, TOPAZ |2|] and AMY 
0. Due to the increased cm. energy of the TRISTAN ring as compared to earlier low energy 
experiments at PETRA and PEP the TOPAZ and AMY cross sections extend to transverse 
momenta of up to 8 GeV. It is hoped that coming data from 77 collisions at LEP1 and 
particularly from LEP2 will enlarge the px range further. Production of high-pr jets probes 
the short-distance dynamics of photon-photon reactions. In addition to providing tests of 
perturbative QCD, data from 77 reactions give us information on the photon structure function. 

In leading order QCD (LO) three distinct classes of contributions to the cross section |4[] 
are identified: (i) The direct contribution, in which the two photons couple directly to quarks 
and neither of the photons is resolved into its partonic constituents (DD component in the 
following), (ii) The single-resolved contribution, where one of the photons interacts with the 
partonic constituents of the other (DR component), (iii) The double-resolved contribution, 
where both photons are resolved into partonic constituents before the hard scattering process 
takes place (RR component). In the DD component we have only the two high-p^ jets in the 
final state and no additional spectator jets. In the DR case one spectator jet coming from low 
transverse momentum fragments of one of the photons is present and in the RR component we 
have two such spectator or photon remnant jets. Experimental separation of the three classes 
is in principle possible and depends whether the experimental arrangements allow to tag the 
spectator jets. 

In next-to-leading order QCD (NLO) the photon-quark collinear singularity arising in the 
DD and DR components is subtracted and absorbed into the photon structure function in 
accord with the factorization theorem. This subtraction procedure at the factorization scale M 
introduces an interdependence of the three components so that a unique separation into DD, 
DR and RR classes is not possible anymore. This means that in NLO all three components 
must be considered together and consistently be calculated up to NLO in the photon structure 
functions || || and in the hard scattering cross sections. 

The photon structure functions are inherently non-perturbative quantities whose magnitude 
and dependence on the fractional momentum x of the outgoing parton must be measured at 
a reference scale M 2 . The change with M 2 is obtained from perturbative QCD evolution 
equations. 

Complete NLO calculations for high-pr jet production in 77 reactions have been done 
previously for the inclusive single jet cross section |], [7| and compared to experimental data 
from TRISTAN 0, ||. The double resolved contribution to the NLO single jet inclusive cross 
section has been investigated in || for TRISTAN, LEP1 and LEP2 energies in order to study 
the dependence on the photon structure function. The restriction to the RR contribution was 
motivated in || that by tagging the two spectator jets it could be isolated experimentally from 
the other components. Unfortunately this is possible only in LO since NLO corrections to the 
DD and DR components can produce spectator jets due to the separation of singularities with 
the help of kinematic constraints. 

In our previous work Q we also calculated the complete NLO inclusive dijet cross sections 
for the DD and DR contributions and estimated the RR contribution by a LO calculation with 
k factors taken from the NLO inclusive single-jet cross section. This was justified since at 
TRISTAN energies, for which the calculation was done, the RR cross section contributes only a 
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small fraction of the total sum in the high-p^ region. In this paper we extend these calculations 
to LEP1, LEP2 and NLC collider energies. Although at these higher energies and at moderate 
Pt the double resolved cross section is much more important, we expect that the estimate of 
the RR dijet cross sections with a k factor will still be reliable. In addition we shall present 
complete NLO predictions for the single-jet inclusive cross section at these energies which we 
expect to be measured first. Since in our earlier work where we compared also to experimental 
dijet cross sections measured by TOPAZ and AMY no details for the calculation of the cross 
sections were given, we shall present them in this longer paper. 

The outline of the rest of the paper is as follows. In section 2 we first give the DD cross 
section in leading order also to fix the notation. Then we present the details of the next-to- 
leading order calculation for the DD two-jet cross section in section 3. 

The calculation of the DR and RR cross sections is based on the work of |T(| and [|TTJ . In 



10|| the formalism for inclusive one- and two-jet cross sections is worked out and applied to the 



calculation of the direct process on jet production in low Q 2 ep collisions. In ||TTJ the inclusive 
one-jet cross section for ep collisions with a resolved photon was calculated which can easily 
be applied to the RR cross section in 77 reactions. 

The numerical results for the DD, DR and RR one- and two-jet cross sections are given in 
section 4. In section 5 we summarize our results and draw some conclusions. 



2 Leading Order Cross Section 
2.1 Photon Spectrum 

The cross sections which we have computed are for kinematical conditions as we expect them 
for 77 collisions at LEP1, LEP2 and the NLC. For LEP1 and LEP2 the spectrum of the virtual 
photons is described in the Weizsacker-Williams approximation (WWA) by the formula |12| 



„ , , a l + (l-x a )\ E 2 ei(l-x a y + m 2 e x 

2lt [ X a m e X a 



+2(1 -x a ) 



mix a 1 



E 2 ®l(l - Xa) 2 + m 2 e x\ x a 



(1) 



where m e is the electron mass and x a = {q a k b ) / (k b k a ) ~ E^/E is the fraction of the initial 
positron energy transferred to the photon with E being the beam energy of the incoming 
positrons and electrons, respectively. k a (k b ) is the four momentum of the incoming positron 
(electron) and q a = k a — k' a is the four momentum of the virtual photon emitted by the positron. 
The momentum fraction on the electron side is denoted as Xb with x b = (k a q b ) / \k a k b ) and 
q b = k b - k' b . 

In the equivalent photon approximation, the cross section for e + + e~ —>■ e + + e" + X with 
arbitrary final state X is then given by the convolution 



1 1 



do- (e + + e -»• e + + e + X) = / dx a J dx b F la / e (x a )F lb / e (x b ) dcr(7 a 7 6 -> X), (2) 



^■a,min ^&,min 



where da{pf a ^ b — > X) denotes the cross section for 7 a 7 b — > X with real photons of energies 
E la = x a E and E Jb = x b E, respectively. 
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The LEP1 cross sections are calculated with v^S* = 2E = 90GeV. The angle G c is the 
maximum angle under which the positrons (electrons) are tagged, which we assume G c = 3.0°. 
For the LEP2 cross section we have chosen v^S* = 175 GeV and C = 30 mrad. For the NLC 
collider option the photon spectra will be specified later when we present our results. 



2.2 Jet Cross Sections in LO 

In this section we write down the leading order cross section for the production of jets with two 
direct photons, i.e. the DD cross section in LO. This is given by the differential cross section 
for two-jet production in the quark parton model, i.e. the cross section for 'ja'jb — > QQ, which 
must be convoluted with the photon distribution functions according to (0). 

The final state quark (antiquark) has momentum pifa) which can be expressed by their 
transverse momentum px and rapidities 771 and r\i- The convention is that the z direction is 
parallel to the electron beam direction. From energy and momentum conservation one obtains 

*a = f§ (e"* + e^) (3) 
x b = H^+e^ 2 ). (4) 

Thus, the kinematical variables of the two jets are related to the scaling variables x a and Xb- 
Under the actual experimental conditions x a and Xb are restricted to fixed intervals, x m i n < 
x a ,Xb < x max < 1. We shall disregard these constraints and allow x a and Xb to vary in the 
kinematically allowed range x at b,mm < x a ,Xb < 1, where 

_ p T e~^ 

a ' mm 2E-p T e*n {> 

and a similar equation for Xb, m m- From eq. fl3|) and (f|), we can express Xb as a function of Pt, 
r\\ and x a : 

x a PTe Vl 

x b = F - „- Vl - ( 6 ) 



2x a E - p T e 



Depending on x a , different regions of the photon energy Xb contribute. Whereas for fixed x a 
and Xb, i.e. fixed energies for both photons, rji and rj 2 are fully determined by px (up to sign 
ambiguities), 771 and 7/2 are allowed to vary due to the kinematical range of the photon energies 
in the intervals x a ^ € [x a ^ m im !]• ^a,6,min are obtained from (5) and are equal to px 2 / ' E 2 . 
The two-jet cross section for e + + e~ — > e + + e~ + jet x + jet 2 is obtained from 

dprdTncfya = XaF ^/^ XbF ^/^ Xb ^ ^ b -^PiPa)- ( 7 ) 

— stands for the differential cross section of the process 'jalb V\V%- The invariants of this 
at 

process are s = (q a + qb) 2 , t — (qb — p\) 2 and u = (qb — P2) 2 ■ They are expressed by the final 
state variables pr, rji and r\i and the initial state momentum fractions x a and Xb- 

s = 4x a x b E 2 (8) 
t = -2x a Ep T e m = -2xbEp T e- Vl (9) 
u = -2x a Ep T e Vl = -2x b Ep T e- m . (10) 
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So, the dependence of the two-jet cross section on p T , rji and r) 2 is determined through the two 
photon distribution functions and the cross section for the 77 subprocess, which depends on s, 
t and u. 

For the inclusive one-jet cross section, we must integrate over one of the rapidities in (0). 
We integrate over 772 and transform to the variable x a using (|^). The result is the cross section 
for e + + e~ — > e + + e~ + jet + X, which depends on p T and 77: 

dV /" 4£'p r dcr 

d^ = 7 2XqE _ pre -, ^ (7a7> - P1P2) ■ (11) 

min 

Here, is given by @ with r/i = 77. 

The cross section for the subprocess 7 a 7b — > q^q\ is well known and is given by 

dt = — ^(i + i) (12) 

The index i denotes the quark flavour and the quark charge. Nq is the number of colours. 
All quarks are considered massless. 



3 Next— To— Leading Order Cross Sections 

The next-to-leading order corrections are calculated with the help of dimensional regulariza- 
tion. In the following subsections we shall consider the virtual and the real corrections needed 
to obtain a finite cross section for the DD case in the limit n — > 4. 



3.1 Virtual Corrections up to 0(a 2 a s ) 

The one loop diagrams for 77 — > qq have an additional virtual gluon, which leads to an extra 
factor a s . These diagrams must be multiplied with the LO diagrams to produce the virtual 
corrections to the 2 — > 2 cross section up to 0(a 2 a s ). These corrections are well known for many 
years now |TB[ . The result can also be obtained from the virtual correction for 75 — > gq [[TTJ 
by taking only the diagrams without the three gluon vertex and substituting the appropriate 
colour factors. The result can be written as 

Mrt-+wk) = eW £ (^) ^ ra - 2 e ) 2NcCFV ~< {s ^ u) + Q(e) ' (13) 

Hy gives the virtual correction to the corresponding reactions up to the n-dimensional phase 
space factor 

dPS (2) 1 /4tts\ . 

(14) 



dt T(l - e) V ut J 8ws 

and the flux factor l/(2s). The expression for V 7 (s,t,u) can be found in appendix A. The 
singular terms oc 1/e 2 and 1/e are proportional to the LO cross section 
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t U . 

- + - (1-e -2e 
u t 



(15) 



with e = (4 - n)/2. 
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3.2 Real Corrections up to 0(a 2 a s ) 

For the 2^3 contributions up to 0(a 2 a s ), we have to take into account all diagrams with an 
additional gluon in the final state. The diagrams are shown in Fig. [TJ 

The four-vectors of these subprocesses will be labeled by q a qb — > PiP2P3- The invariants 
will be denoted by Sjj = (pi + Pj) 2 , = a, b, 1,2,3). For massless partons, the 2^3 
contributions contain singularities at = 0. They can be extracted with the dimensional 
regularization method and can be cancelled against those which originate from the one-loop 
contributions or are absorbed in the renormalized photon structure functions. 



To achieve this, we go through the same steps as described for example in [10). First, 
we calculated the 2 — > 3 subprocesses in n dimensions. In our case we have two classes of 
singularities. Examples of these are shown in Fig. ||. The X marks the propagator leading 
to the divergence. In the first graph, the vanishing of the invariant S12 leads to a final state 
singularity and the second graph becomes singular for s& 3 = 0, which then leads to an initial 
state singularity. 

When squaring the sum of all diagrams in Fig. [| we encounter terms where more than one 
of the invariants become singular, e.g. when the gluon momentum pi — > 0, so that = and 
S13 = 0. These infrared singularities are disentangled by a partial fractioning decomposition, 
so that every term has only one vanishing denominator. Non-singular terms are discarded. It 
turns out that in this limit the results are always proportional to the LO cross sections involved 
in the hard scattering, namely T 7 and T q , where T q stands for 'jq — > gq, so that they can be 
written as 

H FJ = KpjT^. (16) 

Here F and I denote the contributions originating from the final(F) and initial(I) state singu- 
larities, respectively. 

The last step is to integrate the decomposed matrix elements analytically in the region 
Sij < ys, where y is the upper bound for the integration and has to be chosen small enough so 
that it is justified to neglect terms of 0(y), which was assumed already in the partial fractioning 
above. This means we use an invariant mass cut to separate the genuine 2^3 contributions 
from the 2-^2 contributions. For the terms with final state singularities y determines the 
phase space region where two partons, i.e. a quark and a gluon, are combined in one jet. In 
the case of initial state singularities y determines the boundary between the remnant jet of 
the photon and the 2 — > 2 hard scattering 7g — > gq. In || the separation of the final state 
singularities was done with the Sterman-Weinberg parameters e and 5, where e is the cut on 
the gluon energy and S measures the angle between the momenta of the quark and the gluon in 
the recombination of a quark and a gluon into one jet. The integration produces terms oc 1/e 2 
and 1/e, which cancel against those in the virtual corrections or be absorbed into the photon 
structure functions. In the following we shall give the results for the final and the initial state 
singularities separately. 

3.2.1 Final State Singularities 

In this subsection, we assume that after partial fractioning the 2 — > 3 matrix elements are 
singular only for s 12 = 0. For the integration over the singular phase space region we choose 
as coordinate system the cm. system of the partons p\ and p2- The angles of the other parton 
three-momenta q a and p% with respect to p\ and P2 are shown in Fig. [| x is the angle between 
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the momenta q a and p 3 , 9 is the angle between q a and p±, and is the azimuthal angle between 
the planes defined by q a and p\ and p a and p 3 , respectively. Instead of 9 we also use the variable 

b = ^(l-cosfl). (17) 

The angle can be integrated out easily, because the matrix elements do not depend on it in 
the limit that non-singular terms are discarded. 
Here, we define the invariants 

s = (q a + q b ) 2 (18) 
t = (q b ~ Pi ~ P2) 2 ~ 2 Pl p 2 (19) 
u = (qt-ps) 2 (20) 

which differ from the corresponding two-body invariants, but are equal to them for p\ = or 
p 2 collinear with p ± . The variable to be integrated is 

z> = ™. (21) 
q a -q b 

The three-body phase space in n dimensions can be factorized into 

dPS (3) = dPS (2) dPS (r) , (22) 



where 



and 



with 



dPS (2) 1 /Atts 



dt T(l - e) V tu J 8tts 



(23) 



dUF = dzV-h+^j ^(l-6)-^sin- 2 ^. (25) 



Nb and N$ are normalization factors: 

1 





7T 



N t = y d ^sm- 2 ^ = 4%SLJl|. (27) 

^ ' 

The full range of integration is given by z' G [0, — t/s], b G [0, 1] and (j) G [0, 7r]. The singular 
region is defined by the requirement that partons p\ and p 2 are recombined, which means 
S12 — > 0. We integrate over this region analytically up to s 12 < ys, which restricts the range of 
integration to < z' < min{— t/s,y} = yp. 

In the 77 case we obtain the final state matrix element in the following form 

'6PS^H F = Qt^ € (—) ~o I WcCfF^s, t, u) + O(e) . (28) 



27iT(l-2e)' 

The factor F 7 (s, t, -u) can be found in Appendix B. It contains infrared and collinear singularities, 
which cancel against those in the virtual corrections. 



3.2.2 Initial State Singularities 

Here we integrate over the singularity s b s = 0, we use the same p\-p2 cm. system as in the last 
section. In this case, the outgoing quark p^ is collinear to the incoming photon momentum q b . 
So it becomes part of the photon remnant and p\ and P2 are the momenta of the final state in 
the 2-body subprocess. We introduce the new variable 

z b = ™*e[X b ,l], (29) 
q a qb 

where X b = (P1P2) / (lakb) — E q j 'E is the fraction of the initial electron energy transferred to 
the quark. 

We have now the following definition for the Mandelstam variables 

s = (q a + z b q b ) 2 (30) 
* = (qa-pi) 2 (31) 
u = {q a -P2? (32) 

and the variable, which parametrizes the singular region, is now 

z" = (33) 

q a qb 

Again the three-body phase space factorizes into 

dPS (3) = dPS (2) dPS (r) , (34) 

where dPS^ 2 -* is again the phase space of the 2^2 process now in the limit q b p^ — > given by 
(|23|), and 

dps< " - (-^W^Tt^''^' < 35 > 



s 



where 



„ „^ e dz b ( z b y d<f) . _ 2e / l - 26 

dm = dz z — — sin 0— 36 

z b \l-z b J T 2 (l-e) 



and 



-l+2e 

H a (z") = + (1-T^— ) =± + 0(z") (37) 



z \ I z 



can be approximated by 1, because it leads only to negligible terms of 0(y). The full region of 
integration is given by z" G [0, —u/s], z b G [X b , 1] and (p G [0,7r], where again the dependence 
on (f) can easily be integrated out. The singular region where the integration is done analytically 
is given by the requirement < z" < min{— u/s, y} = yj. 
The result is 

/dPSMff, = / ^gf/x 4e (^) £ £^Y^Y^ a ^A(^; S ,t, W ) + 0(e). (38) 



A',, 
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Again, I y (z; s,t,u) can be found in the appendix. 

In the case of initial state singularities the following quark-photon matrix element factorizes 



TW(M,«) = (l- £ ) 



s t. 



The case z" = q a Vzl lalb leads to the same result with (z b 
explicitly the pole in 1/e proportional to 



N c Ql 2z 2 -2z+l 



2e . (39) 
-> z a ). The term J 7 shows 

(40) 



This function appears in the evolution equation of the photon structure function as an inho- 
mogeneous or so-called point-like term. Therefore, the photon initial state singularities can be 
absorbed into the photon structure function. 



3.3 DD Jet Cross Section in NLO 

To obtain a finite cross section for the DD case, we must add the parts considered in section 
3.1 and 3.2. Then the poles in 1/e and 1/e 2 cancel and we can take the limit e — ► 0. The result 
is a special kind of two-jet cross section, where the recombination of two partons into one jet 
or the recombination of a parton with the photon remnant jet is done with an invariant mass 
cut-off y. Including the NLO corrections, we get 



dV 



dpTd77idr72 



■Ea Fy a /e (•£ a ) bF'yi, /e ( b ) 



da 



(7a7b -> pm) 



(41) 



In (0), -^(7a7& P1P2) stands for the two-body contribution in LO and NLO together 
with analytically integrated contributions of the soft and collinear divergent regions of the 
three-parton final state. The contributions from the initial state singularities are denoted 

~r (laQ P1P2) and —{q^b —> P1P2), where either the photon 7& or 7 a has collinear singular 

Lit- Lit- 

terms removed, respectively. 



The two-body contribution can be written as 
da 



dt 



{lalb -> PlP2) 



2tt 



-C(T 7 A 7 + B. 



1)1 



(42) 



where 



B- 



T 



C F 
C F 



-7T 2 + In 2 — + In 2 — - 2 In 2 y F - 3 \ny F 

.3 s s 

—t —u u —t t —u 

2 In — + 2 In — + 3- In — +3- In — 

s s t s u s 



u 



-u 



2 + - ln z + 2 + - ln z — 



t 



u 



-t 



- + - ) , and C 

t u 



Qt N c- 



(43) 



(44) 
(45) 
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The cut dependence of the final state corrections in A 7 is contained in the yp dependent 
terms. For yp — > 0, these terms behave like (— In 2 yp), which leads to unphysical negative cross 
sections for very small yp. Thus if y is used as a physical cut, it must be sufficiently large. In 
most applications, we shall use these results for computing inclusive cross sections, in which the 
y dependence of the two-jet cross section cancels against the y dependence of the numerically 
calculated three-jet cross section. 

The two-jet cross section for the initial state is calculated from 



da 
dt 



{lib 



da 
~dt 



PlP2j 



(lag -> PlP2) 



da a 



dZg 

Za 



Pq^j(Za) In 



X a 

da a 
-{laq^pm) 



dt 

1 



2n 



dz h 



Pq^j(zb) In 




(1 - z a )y!S S 



1 - z h )yis 



z h Ml 



-l)+N c Q 2 



NcQl 



(46) 



(47) 



In (|H) and fl47f ), M a resp. Mb are the factorization scales. The dependence on M a resp. Mb 

must cancel against the M a resp. Mb dependences of the LO DR and RD contributions. The 
da 

cross section — has the following form 
dt 



da 

Tt {iq 



Q9) 



'-Q-C, 



(48) 



All the results of this section are for the MS subtraction and renormalization scheme. 



4 Inclusive One— and Two— Jet Cross Sections 

In this section we present some characteristic numerical results for one- and two-jet inclusive 
cross sections which have been obtained with our method of slicing the phase space with in- 
variant mass cuts. In a short communication we used this method to calculate the differential 
one- and two-jet cross section as a function of px integrated over special rapidity intervals 
and compared it to recent experimental data of the TOPAZ || and AMY |J collaborations 
at TRISTAN. For this paper we have calculated various one- and two-jet distributions with- 
out applying special cuts on kinematical variables of the initial or final state dictated by the 
experimental analysis, although our approach is particularly suitable for this. 

The calculation of the cross sections proceeds as follows. For the DD and DR components 
we use the phase space slicing method for the inclusive one- and two-jet cross sections. The 
details for the calculation of the DD component below the cut are described in the previous 
section. The DR component is identical with the work in [10| for the direct photoproduction 



of jets in low Q 2 ep collisions. We must replace the proton structure function by the photon 
structure function and obtain the cross sections for the single-resolved contribution. With 
the results presented in [IU| and in the previous section we are able to calculate the inclusive 



cross section for one- and two-jet production. For the double-resolved contribution the NLO 
corrections for the two-jet cross section are not available yet. Here only the LO cross sections 
are at our disposal. To estimate the NLO corrections we make use of the NLO cross sections 
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for the inclusive one-jet case. This cross section has the same structure as the NLO resolved 
inclusive one-jet cross section for photoproduction — > jet + X, where the photon structure 
function at one vertex is replaced by the proton structure function. For this cross section the 



NLO corrections are known from earlier work |[L1| . We use these results and transform the 
cross section to the 7 case by replacing the proton structure function by the photon structure 
function. This gives us the full NLO inclusive one-jet cross section for which we present results 
later. We compare the NLO cross section with the RR LO cross section and calculate the k 
factor from it. Then the same k factor is applied to the LO two-jet cross section. We checked 
this procedure with the DD and DR cross sections, where we know also the NLO two-jet cross 
sections. For these cases we found that the k factors for the one- and two-jet cross sections 
are approximately equal. It is reasonable to expect this also to be the case for the RR cross 
sections. 

The further calculation of the DD and DR cross sections is based on two separate contri- 
butions — a set of two-body contributions and a set of three-body contributions. Each set is 
completely finite, as all singularities have been cancelled or absorbed into structure functions. 
Each part depends separately on the cut-off y. If y is chosen large enough, the two parts 
determine physically well defined two-jet and three-jet cross sections. However, our analytic 
calculations are valid only for very small y, since terms 0(y) have been neglected in the analytic 
integrations. For very small y, the two cross sections have no physical meaning. In this case, 
the (In?/) terms force the two-body contributions to become negative, whereas the three-body 
cross sections are large and positive. When both contributions are added to yield a suitable 
inclusive cross section, as for example the inclusive one-jet cross section, the dependence on 
the cut-off y will cancel. Then, the separation of the two y dependent contributions is only a 
technical device. The cut-off only serves to distinguish the phase space regions, where the inte- 
grations are done analytically, from those where they are performed numerically. Furthermore, 
y must be chosen sufficiently small so that experimental cuts imposed on kinematical variables 
of the final state do not interfere with the cancellation of the y dependence. 

First we consider the inclusive one-jet cross section. We choose the definition of the Snow- 
mass meeting |l4j] for combining two nearly collinear partons. According to this definition, 



two partons i and j are recombined if R i: j < R where Ri = y(^j — t]j) 2 + (0j — <pj) 2 is the 
distance between parton i and jet J in the rapidity-azimuthal space. 77^, 0j and t]j, <pj are the 
rapidities and the azimuthal angles of parton i and of the recombined jet J, respectively. We 
choose R = 1 in all the following results. The Snowmass condition means that two partons are 
considered as two separate jets or as a single jet depending whether they lie outside or inside 
the cone with radius R around the jet momentum. Unfortunately this definition is not unique. 
In some cases it may happen that two partons i and j qualify both as two individual jets i 



and j and as a recombined jet ij. In this case we count only the combined jet following |L5 
In NLO the final state may consist of two or three jets. The three-jet sample consists of all 
three-body contributions, which do not fulfill the cone condition. The inclusive cross section 
will depend on the value of R chosen. It will increase with increasing R. 

Before we calculate the final results to be presented in the figures we have made some checks 
of the NLO corrections to the one- and two-jet cross sections. First we checked that the DD 
and DR cross section are independent of the slicing cut y if y is chosen small enough. This 
was the case for y < 10~ 3 in all considered cases. For y > 10~ 3 we observed some small y 
dependence which is caused by our approximation that we neglected contributions 0(y) in the 
analytical contributions to the two-jet cross section. In the final evaluation we use y = 10~ 3 . 
Furthermore we tested that the sum of the NLO direct and the LO single resolved cross section 
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is independent of the factorization scale M. The same test was performed for the sum of the 
NLO single resolved and the LO double resolved cross section. These tests were done for the 
one- and two-jet cross sections separately. Similar checks have been reported for the one-jet 
photoproduction cross section flff] . 

The input for our calculation is as follows. For the DR and RR cross sections we need the 
parton distributions, Fj/ 7 , in the photon. We have chosen the NLO set of Gliick, Reya and 
Vogt (GRV) in the MS scheme 0. This means that the DD, DR and RR cross sections must 
also be calculated with the MS subtraction. We choose all scales /i = M = pt and calculate 
a s (fi) from the two-loop formula with Nf = 5 massless flavours with = O.lSOGeV^ equal 
to the A value of the NLO GRV photon structure function. The charm and bottom quarks are 
treated also as light flavours with the boundary condition that the charm (bottom) contend 
of the photon vanishes for M 2 < m 2 (ml) (m c = 1.5 GeV, m^ = 5 GeV). We do not apply a 
special cut on the energy fractions x a and Xb in (H) but integrate from x ajm i n and x^min to 1, 
where x a>m ^ n and x^min are given by kinematics (see (IP). 

In the following we show results for the three different cm. energies: (i) LEP1 with = 
90 GeV and 6 C = 3° in the photon spectrum in <g), (ii) LEP2 with yfs = 175 GeV and 
6 C = 30 mrad and (hi) NLC in the TESLA design with \/~S = 500 GeV and the photon spectra 
given by the sum of the WWA spectrum in (|I|) with 9 C = 175 mrad and the beamstrahlung 



spectrum given in with parameters Y e g = 0.039 and o z = 0.5mm []18 |. (See |i| for the 
impact of various collider options.) 

First we show the inclusive one-jet cross sections d 2 a/dpxdr] for the three machines. In Fig. f| 
this cross section is plotted as a function of pr for rapidity 77 = 0, where this cross section is 
maximal. Only the NLO predictions are plotted for the DD, DR and RR cross sections and for 
the sum of all contributions. Below pt = 5 GeV the RR component is dominant whereas for 
the larger p? the DD cross section gives the largest contribution. At p? = 5 GeV the DD and 
RR cross sections are equal. Above pt = 10 GeV the DR cross section is larger than the RR 
cross section. The rapidity distribution at px = 5 GeV is shown in Fig. |5|, again for the three 
components and for the sum. The slight variations in the RR curve are caused by the limited 
accuracy of the numerical integrations. All three cross sections must be symmetric at 77 = 0, 
where they exhibit a rather broad plateau. The same distributions, i.e. d 2 a/dpTdr] for 77 = 
and d 2 a /dpTdr] for px = 10 GeV are shown in Fig. || and [7] for LEP2. The qualitative behaviour 
of the one-jet cross sections is similar. Due to the larger cm. energy the cross sections for LEP2 
are larger than for LEP1. The RR contribution is more significant now. The DD distribution 
crosses the RR distribution at larger p? than in the LEP1 case. The rapidity distribution at 
Pt = 10 GeV looks very similar to that in Fig. |5|. At this pt value all three components make 
significant contributions to the total sum. In the NLC case the pattern is somewhat different 
as can be seen in Fig. ^| and |9]. Due to the further increase of the cm. energy the cross section 
is increased by an order of magnitude as compared to the LEP2 cross section. The hierarchy 
between the DD, DR and RR components has changed somewhat. The RR cross sections lie 
for pt > 5 GeV below the DD and DR cross sections. All three components cross each other 
below pt = 5 GeV, where RR starts to dominate. Above pt = 5 GeV the pt distribution is 
dominated by the DD component. The total rapidity distribution at pt = lOGeV^ is less flat 
around 77 = 0. This originates from the DD and DR part. The DD part has the steepest 
behaviour towards the kinematic boundaries. The different functional behaviour of da/dpT 
and d<r/d77 in Fig. [5] and [] as compared to the cross sections for LEP2 comes from the changed 
photon spectra which is now a superposition of the WWA and the beamstrahlung spectrum. 
The rapidity distribution for the NLC looks similar to the results obtained recently for charm 
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quark production in two-photon collisions [30]. However, the difference between the LEP2 and 
NLC rapidity distributions is even more drastic for heavy quark production than in Fig. |7| and 
B[ By comparing with the LO one-jet cross section, where the same photon structure function 
and the same a s as in the NLO cross section is used we obtained the k factors for the RR 
contributions. They are k = 1.85, 1.90 and 1.90 for the three cases LEP1, LEP2 and NLC. We 
shall use these k factors for correcting the LO RR predictions of the two-jet cross sections. 

Next we consider the results on inclusive two-jet production which will be shown in the 
following figures for the LEP2 case. In Fig. [U], [11] and [12|, we present dV/dpyd^idr^ as a 
function of pr l for rji = and various choices of 772 = 0, 1, 2. Here, pt x and rji are the transverse 
momentum and the rapidity of the so-called trigger jet. 772 is the rapidity of the second jet, 
so that an d Pt 2 are the two highest transverse momenta of the three-jet configuration. For 



exactly two jets in the final state, we have pt x = Pt 2 ■ in Fig- O an d we can see how the cross 
section decreases when 772 is chosen away from the maximum region at 772 = 0. In particular the 
RR component becomes more important for increasing pt when 7/2 increases, 771 — is always 
kept fixed. Comparing the cross sections in Fig. |l^ and |TT] we see that the sum of the three 
components hardly changes when 772 = 1 instead of 7/2 = 0, although the relation of the DD, 
DR and RR cross sections is different in the two cases. We have studied the inclusive two-jet 
cross section also as a function of rjx and r/ 2 for fixed pn- As an example, we show for LEP2 
the two-dimensional distribution dV/dprdr/idr^ for p^ = lQGeV in form of a lego-plot in 
the intervals 771, 772 G [—3.0, 3.0]. The cross sections for the DD, DR and RR contributions and 
for the sum are shown separately. The DR cross section contains both contributions with the 
resolved photon at the upper or the lower vertex. If they are considered separately we checked 
that they are symmetric for r/i <-> rj 2 . 

For the NLC case we present only dV/dprdr/idr^ as a function of for r/i = 7/ 2 = 
where the two-jet cross section is maximal. This is shown in Fig. O. Comparing with the 



corresponding cross section for LEP2 in Fig. 10 we notice the increase of the cross section 



by more than a factor 10 and that the RR component is somewhat more significant also for 
larger p T . The cross section dV/dp^dr/xd?^ as a function of t/ 2 for p Tl = 10 GeV and r)i — is 
presented in Fig. pj^ for the three components DD, DR, RR and the sum. Only the RR cross 
section is rather flat near the maximum 772 = 0. The other components are much steeper away 
from 772 = as it is the case also in the sum similar to the one-jet cross section shown in Fig. || 
It is clear that many more distributions or partially integrated cross sections using other 
two-jet variables can be calculated with the phase space slicing method. So, for example, one 
could study kinematic regions where either the DD or the RR cross section is enhanced, as has 
been done for two-jet production in 7p processes JlpJ by making cuts in x 7 , the fraction of the 
photon energy participating in the hard scattering process. Other interesting topics are the 
cone dependence of the inclusive cross sections or the study of the invariant mass distribution 
of the two jets for different rapidities or the angular distribution of the two jets to test the 
parton-parton scattering dynamics in a different way. 



5 Summary and Conclusions 

Various inclusive one- and two-jet cross sections have been calculated for the direct, single 
resolved and double resolved contributions as a function of px and jet rapidities in NLO. For 
the double resolved two-jet cross section the NLO corrections are estimated with a k factor 
taken from the inclusive one-jet cross section. For the direct and single resolved components 
all NLO corrections are fully evaluated. Infrared and collinear singularities are cancelled with 
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the phase space slicing method using an invariant mass cut-off. This method is particularly 
useful for incorporating cuts on the final state and for obtaining results with different choices 
of jet algorithms. Analytical formulas for the different contributions giving the dependence 
on the slicing parameter are derived for the direct contribution. The same results for the 
single resolved contribution are taken from the corresponding calculation of jet production in 
7P processes 
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Numerical results for the inclusive one-jet and two-jet cross sections in three energy ranges 
corresponding to LEP1, LEP2 and NLC have been presented. In the NLC case we have calcu- 
lated the photon spectra from a superimposition of bremsstrahlung and beamstrahlung spectra. 
For this case the photon spectra lead to quite different rapidity distributions as compared to 
LEP1 and LEP2 where we have only the bremsstrahlung spectrum. The hierarchy of direct, sin- 
gle and double resolved cross section is very similar for the three machines. The RR component 
dominates only for rather small < 5 GeV in all cases. For larger p? the direct contribution 
is dominant. 

Since our results have been successfully tested already at smaller energies by comparing 
with one- and two-jet measurements of the TOPAZ and AMY collaborations we are confident 
that our results are also reliable for the larger energies achievable at LEP2 and NLC. 
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A Virtual Correction 



Vy(S, t, U) 



2 3 2tt 2 



7 + l n 2^ + ln 2^ 

S 



T<?\s,t,u) 



^ , — t n , —u 
+2 In h 2 In 

s s 



3^1n^ + 3*ln^ 

t s us 



_( 2 + -)ln 2 — + f 2 + -)ln 2 — 



(49) 



B Final State Correction 



F 1 (s, t, u) 



2 3 

e 2 e 



+ 7 — 2 In 2 — 3 In yp 



(50) 



C Initial State Correction 



Ij(z;s,t,u) = 



1 1 



; P q ^(z) + 1 + In (^—^j (2z 2 -2z + l) Tf\s,t,u) (51) 



eN c Q 2 ' 

Here one needs the Altarelli-Parisi splitting function for the process photon — > quark, which is 



P q ^(z) = N C Q 2 



2z 2 + 2z+l 



(52) 
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4 Figure Captions 

Fig. 1: The three body diagrams for 77 — » qqg. 

Fig. 2: Three body diagrams with final and initial state singularities. 

Fig. 3: Kinematic diagram for the three-body final state defining the angles in the cm. system 
of partons p\ and p 2 . 

Fig. 4: NLO inclusive one-jet cross section d 2 a/dprd?7 as a function of px at i] = for direct 
(DD), single resolved (DR) and double resolved (RR) contributions and for the sum of 
all components. LEP1 photon spectra, v^S* = 90GeV. 

Fig. 5: NLO inclusive one-jet cross section as a function of rj at p T = 5 GeV and for DD, DR, 
RR and the sum. LEP1 photon spectra, V~S = 90 GeV. 

Fig. 6: Same as Fig. 4. LEP2 photon spectra, y/S = 175 GeV. 

Fig. 7: Same as Fig. 5, p T = 10 GeV. LEP2 photon spectra, VS = 175 GeV 

Fig. 8: Same as Fig. 4. NLC-TESLA photon spectra. y/S = 500 GeV. 

Fig. 9: Same as Fig. 5, p T = 10 GeV. NLC-TESLA photon spectra. = 500 GeV. 

Fig. 10: NLO inclusive two-jet cross section as a function of pt for rji = 772 = for direct (DD), 
single resolved (DR), double resolved (RR) and sum. LEP2 photon spectra. v^S* = 
175 GeV. RR component estimated with k factor. 

Fig. 11: Same as Fig. 10 for r]i = 0, r] 2 = 1. 

Fig. 12: Same as Fig. 10 for i]i = 0, i] 2 = 2. 

Fig. 13: NLO triple differential cross section dV/dpTdr/idr^ for px = 10 GeV as a function of 771 
and 772 for DD, DR, RR components and the sum. LEP2 photon spectra. =175 GeV . 

Fig. 14: Same as Fig. 10. NLC-TESLA photon spectra. VS = 500 GeV. 

Fig. 15: NLO inclusive two-jet cross section as a function of r] 2 for p T = 10 GeV and % = 0. 
NLC-TESLA photon spectra. y/S = 500 GeV. 
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